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Abstract 

A program for the calculation of the total Higgs production cross section via 
gluon fusion at hadron colliders including next-to-leading order QCD corrections is 
presented. It is suitable especially for Standard Model Higgs bosons and the neutral 
Higgs particles of the minimal supersymmetric extension. The program provides a 
model-independent calculation for scalar [CP-even] as well as pseudoscalar [CV- 
odd] Higgs bosons including the contributions of virtual top and bottom quarks 
inside the loop coupled to the Higgs particles. The relevant input parameters can 
be chosen from an input file. As a special case the minimal supersymmetric exten- 
sion of the Standard Model can be investigated. The corresponding couplings are 
implemented including the leading higher order corrections. 
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1 Introduction 



At the future hadron colhder LHC the Higgs boson will be produced primarily via the 
gluon fusion mechanism for the entire relevant Higgs mass range within the Standard 
Model [SA4] as well as the minimal supersymmetric extension [A^iSiSA^] 0]. Recently 
the next-to-leading order QCD corrections to the production cross sections of scalar [CP- 
even] as well as pseudoscalar [CP-odd] Higgs bosons have been calculated They 
are significant for the theoretical prediction of the cross sections leading to an increase by 
up to a factor of two compared to the lowest order results. 

In this paper the program HIGLU^ for the calculation of the total Higgs production 
cross sections including these next-to-leading order QCD corrections will be presented. 
Various relevant input parameters can be chosen from an input file including a flag spec- 
ifying the model. Possible options are the Standard Model, its minimal supersymmetric 
extension and a general Higgs model by initializing the Higgs Yukawa couplings to the 
heavy quarks appropriately. The program includes the contribution of the top and bottom 
quarks in the loop that generates the Higgs couplings to gluons. Within the Standard 
Model as well as in most of the parameter space of the AiSSAi these contributions pro- 
vide an excellent approximation for all cases in practice. Moreover, the program allows 
to calculate the decay widths of Higgs bosons into gluons including next-to-leading order 
QCD corrections. The gluonic decay mode plays a significant role in the intermediate 
mass range at future e~^e~ colliders |Q. 

The source code of the program is written in FORTRAN. It has been tested on a 
VAX-station using the operating system VMS, on different workstations running under 
UNIX and IBM computers with the operating systems CMS and TSO. The numerical 
integration is performed by using the VEGAS-package [Q for integrals of dimension up 
to three. Parton distributions can be attached to the program in any desirable way 
by adjusting the corresponding subroutine as described in section ^. As the standard 
parametrization the program contains the GRV sets [^. 



2 Results 



2.1 pp^ ^ + X 

The hadron cross section of Higgs bosonQ production via gluon fusion gg —>■ ^ = Ti, A) 
including next-to-leading order QCD corrections, can be cast into the form 

a{pp ^<^ + X) = ato + Aal, + Aa* + Aa% + Aa% (1) 

with the lowest order cross sections 

a{pp^^ + X) = a^T^^ (2) 



■^Comments or suggestions are welcome and should be sent to spira@desy.de. 

"^The scalar [CT'-even] Higgs particles will generically be denoted by H, the pseudoscalar [CT'-odd] by 



A and all the neutral Higgs bosons by <i>. 
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The coefficients are 



TC 



128^2: 



Q 



(3) 



They include the Yukawa couphngs Qq normahzed to the SM. couphngs, and the quark 
amphtudes 



(4) 



The function /(r) is defined as 



arcsm 



fir) = 



and the scahng variables are 



log j== - ITT 

1 — yl — r 



T > 1 



r < 1 



(5) 



Arrir 



and 



(6) 



The parameter uiq denotes the heavy quark mass, m<j, the Higgs boson mass, s the total 
center of mass energy squared, Gp the Fermi constant and as the QCD coupling constant. 
The term Ao"*^^ parametrizes the infrared regularized virtual two-loop corrections and 
the terms Act* (z, j = g, q, q) the individual coUinear regularized real one-loop corrections 
corresponding to the subprocesses 

gg ^ ^g, gq ^q, qq ^g 

The expressions for Ao"*^^ and Aa* can be found in Refs. |^ Q . The gluon luminosity is 
defined by 

^ = j'%{x,Q')g{T/x,Q') (7) 

where g{x,Q'^) denotes the gluon density. The natural values to be chosen for the renor- 
malization scale /i of the strong coupling ^^(/i^) and the factorization scale Q of the parton 
densities is given by the Higgs mass m$. 

The program HIGLU calculates the five terms in eq.(|I]) contributing to the total cross 
section separately as well as their sum for all kinds of neutral Higgs bosons $. The 
computation of differential cross sections can be found in Refs. and is not implemented 
in the program. 
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3 gg 



The decay widths of Higgs bosons $ into gluons up to next-to-leading order are given by 



^{^ 99(9), 9Qq) 

with the leading-order expressions 

^LoO< 99) = 

TioiA -> gg) = 



1 + E, 



TC 



911 



(8) 



36v^7r3 
16v^7r3 



T.9^A^{rQ) 
Q 

T.9^A^{rQ) 



(9) 



The coefficient -E*^^ denotes the infrared reg- 



ularized virtual two-loop corrections, E^^g and E^^^ the coUinear regularized real one-loop 



The amplitudes Aq{tq) are defined in eq, 

and -E* 

corrections. The analytical formulae of these contributions can be found in Ref. [^. The 
parameter Np fixes the number of light external flavors produced in the decay $ — gqq, 
which is defined to be equal to the number of flavors contributing to the QCD jS function. 
This definition maps large logarithms into the running strong coupling agifi"^)- The nat- 
ural renormalization scale fi of the strong coupling Q;s(yU^) is given by the corresponding 
Higgs boson mass m$. 



4 Input Parameters 

In addition to the source code of the program HIGLU an input file defined as unit 98 is 
needed, from which the program reads the input parameters. The name of this input file 
can be defined in the ffist OPEN statement of HIGLU. It should be noted that the input 
numbers must not start before the equality signs in each corresponding line. The input 
file contains the following parameters: 

process: integer 

choose the process to be calculated: 
0: gluon fusion gg ^ ^ 
1: gluonic decay ^ gg 

collider: integer 

choose the hadron collider mode for the gluon fusion process: 
0: pp 
1: pp 

This flag is only relevant for the gluon fusion process. 
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energy: double precision 

center of mass energy [in TeV] of the hadron collider for the gluon fusion process 

model: integer 

choose the model, in which the process should be calculated: 
0: SM 

1: M.SS A4 including the dominant two-loop corrections to the Higgs masses and 
couplings 

2: M.SSM. including only the the leading one-loop corrections, which increase 

as the fourth power of the top mass. 
3: any model 

All stop mixing parameters of the A4SSA4 are set equal to zero. 

tanbeta: double precision 

the A4SSA4 parameter tg/3, which is irrelevant if the flag model equals zero or 
three. 

Qb/. double precision 

the Yukawa couplings of the corresponding Higgs boson to the top and bottom 
quarks normalized to the SAi couplings. These couplings are only relevant, if the 
flag for the model is chosen to be three. If the calculation is performed within 
the AiSSAi, these parameters are irrelevant, because the Yukawa couplings are 
calculated from tg/? and the chosen Higgs mass m Higgs- In the SJ^ [flag model = 
0] these couplings are automatically set equal to unity. 

nib/, double precision 

pole masses of the bottom and top quarks [in GeV] . Useful values for these masses 
are given in the sample of the input file in the appendix: rrit — 176 GeV^, rUb — 5 GeV^. 

type: integer 

choose the neutral Higgs boson type: 
1: heavy scalar H 
2: pseudoscalar A 
3: light scalar h 

In the SM. this flag is set equal to unity automatically. If the model is chosen to 
be the M.SSM., all three types of Higgs bosons can be chosen. If the flag for the 
model is set equal to three, only the values 1 and 2 are possible, because the nature 
of the Higgs boson will be characterized by its mass and its Yukawa couplings gb,t 
to the heavy quarks. 

ITT' Higgs'- doublc precision 

the mass [in GeV] of the chosen Higgs boson type. 

loop: integer 

choose one-loop or two-loop formula for the strong coupling in the MS scheme: 
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1: one-loop 
2: two-loop 



choice: integer 

choose the input value of the strong coupling ct^: 
1: as is defined by the value for as{Mz) 
2: as is fixed by the QCD scale 

alphas : double precision 

strong coupling as at scale of the Z boson mass Mz- This parameter is relevant for 
the fiag choice = 1. 

Np: integer 

the number of fiavors, for which the QCD scale A^^"* is given. This parameter can 
be chosen to be 3,4,5 or 6. The QCD scales for different numbers of fiavors are 
computed by using the quark masses rrib^t from the input file and rric = 1.5 GeV via 
the matching conditions of the MS scheme at the thresholds. 

lambda: double precision 

QCD scale A^^ in GeV. 

The parameters A^^ and lambda are relevant for the flag choice = 2. 
n-ext- integer 

number of light external flavors to be included in the gluonic decay width of the 
Higgs boson. If Uext is set equal to 4, the bottom contribution $ — > bbg is subtracted, 
and if next equals three, the part $ — > ccg is removed. This parameter can be chosen 
to be 3,4 or 5. The running strong coupling as{iJ>^) uses the same number of flavors. 

mui^2'- double precision 

parameters defining the renormalization scale // [in GeV] of the strong coupling con- 
stant as in the following way: 

/i = /il rUHiggs + fJ'2 

The variable /Xi is a dimensionless coefficient of the Higgs mass mniggs, and 
denotes a fixed scale [in GeV] . 

QiX- double precision 

parameters fixing the factorization scale Q [in GeV] for the gluon fusion process: 

Q — Qi mniggs + Q2 

Qi denotes a dimensionless coefficient of the Higgs mass niHiggs, whereas Q2 defines 
a fixed scale [in GeV] . 

abserr: double precision 

absolute error for the VEGAS integration. 
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points: integer 

number of points for the VEGAS integration. 10000 points yield a sufficient preci- 
sion, whereas the results using 1000 points with 5 iterations are already reliable at 
the percent level. 

itmax: integer 

maximal number of iterations for the VEGAS integration, 
print: integer 

flag for the print out style of the VEGAS iterations: 

0: no print out 

1: pretty print out 

10: print out as a table 

A sample of the input file is given in the appendix. 

5 Structure Functions 

For the implementation of structure functions the subroutine STRUG has to be changed 
appropriately. This part of the program reads as follows: 

SUBROUTINE STRUG (X,Q, PDF) 

IMPLICIT DOUBLE PRECISION (A-H,0-Z) 

DIMENSION PDF (-6: 6) 

COMMON/FACSC/ I SCHEME 
C. . .X - BJORKEN X 

C...Q - MOMENTUM SCALE (IN GeV) 

C. . .PDF(-6:6) - MATRIX CONTAINING X*P(X,Q) 

C... IPDF = -6 , -5 , -4 , -3 , -2 , -1 ,0 ,1,2,3,4,5,6 
C. . . T_BAR,B_BAR,C_BAR,S_BAR,U_BAR,D_BAR,GL,D,U,S,C,B,T 

C— CHOOSE PROTON STRUCTURE FUNCTIONS AND THEIR FACTORIZATION SCHEME 
C— ISCHEME: 1 

C MSBAR DIS 

ISCHEME=0 

Q2=Q**2 
ISET=2 

CALL PDGRV(ISET,X,Q2,PDF) 

RETURN 
END 
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The line calling the subroutine PDGRV has to be changed, if another parametrization 
should be used. The flag ISCHEME has to be specified according to the factorization 
scheme, which has been used for the parametrization of the parton densities. Note that 
the array PDF has to be generated in accordance with the given convention. 



6 Output 

The output of the program HIGLU is written to a file unit 99, which contains the chosen 
input parameters as well as all results obtained by VEGAS integrations. The VEGAS 
iterations are written to the standard output, if the flag print is set equal to 1 or 10. The 
integrated results contain all individual contributions of the QCD corrections separately 
and their total sum. A sample of the output is presented in the appendix, which can be 
obtained with the input file shown before, if the flag model and the Higgs mass muiggs 
are specifled appropriatly. If the M.SSM. is chosen as the Higgs model, the pseudoscalar 
Higgs mass corresponding to the scalar Higgs masses m-^ is also typed out. 
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7 Appendix 

7.1 Input File 

An example of the input file is given by 

PROCESS: = GG — > H 1 = H — > GG 

PROCESS = 



COLLIDER: = P P 1 = P PBAR 



COLLIDER = 



TOTAL ENERGY: [TEV] 



ENERGY = 14. DO 



MODEL: = SM 1 = MSSM (TWO-LOOP) 2 = MSSM (ONE-LOOP) 3 = ANY 



MODEL = 



TAN (BETA) : (MSSM) 



TANBETA = 1.5D0 



COUPLINGS: G_B = BOTTOM G_T = TOP 
========== (MODEL = 0) 

G_B = l.DO 

G_T = l.DO 



QUARK MASSES: [GEV] 



M_B = 5. DO 

M_T = 176. DO 
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HIGGS TYPE AND MASS [GEV] : 1 = HEAVY SCALAR 2 = PSEUDOSCALAR 
========================== 3 = LIGHT SCALAR 



TYPE = 1 
M.HIGGS = 200. DO 



SCALES : 


[GEV] MU = 
Q = 


MU_1*M. 
Q_1*M. 


.HIGGS 
.HIGGS 


+ MU_2: 
+ Q_2: 


RENORMALIZATION 
FACTORIZATION 


SCALE 
SCALE 


MU_1 
MU_2 

Q_l 
Q_2 


= l.DO 
= O.DO 
= l.DO 
= O.DO 












ORDER OF 


ALPHA.S : 


1 = LO 


2 


= NLO 







LOOP = 2 



DEFINITION OF ALPHA.S: 1 = ALPHA.S (M_Z) 2 = BY LAMBDA (N_F) 



CHOICE = 1 



ALPHA.S (M_Z): 



ALPHA.S = 0.118D0 



LAMBDA.NF: [GEV] (QCD SCALE) 



N_F = 5 

LAMBDA = 0.226D0 



NUMBER OF EXTERNAL LIGHT FLAVORS: (FOR H ~> GG) 



N_EXT = 3 



VEGAS: ABSERR = ABSOLUTE ERROR 

====== POINTS = NUMBER OF CALLS 

ITMAX = NUMBER OF ITERATIONS 
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PRINT = PRINT OPTION FOR INTERMEDIATE VEGAS-OUTPUT 
1 10 

NO OUPUT PRETTYPRINT TABLE 

ABSERR = O.DO 
POINTS = 10000 
ITMAX = 5 
PRINT = 10 



7.2 Output 

7.2.1 Gluon Fusion 
VEGAS : 

ABSERR = . OOOOOOE+00 

POINTS = 10000 ITERATIONS = 5 

GLUON FUSION: GG — > HIGGS 



P P COLLIDER 



ENERGY 



14.0000 



TEV 



LAMBDA_5 = 
REN-SCALE = 



0.226232 GEV NLO-ALPHA_S (M_Z) = 0.118000 

200.000 GEV FAC-SCALE = 200.000 GEV 



T-MASS 



176.000 



GEV 



B-MASS 



5.00000 



GEV 



HIGGS 



= H 



M_H 




200.000 


G"H_B 




1.00000 


SIG.LO 




6 . 42484 


SIG.VIRT 




3.45509 


SIG.GG 




4.78077 


SIG.GQ 




-0.268282 


SIG.QQ 




0.182952E 


SIG.NLO 




14.4107 



GEV 

G"H_T = 1.00000 



+- 


0, 


.681792E- 


-05 


PB 


+- 


0, 


.366648E- 


-05 


PB 


+- 


0, 


.466710E- 


-02 


PB 


+- 


0, 


,709725E- 


-03 


PB 


+- 


0. 


. 173479E- 


-04 


PB 


+- 


0, 


.472080E- 


-02 


PB 
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7.2.2 Gluonic Decay 



VEGAS : 

ABSERR = . OOOOOOE+00 

POINTS = 10000 ITERATIONS = 5 

HIGGS — > GG 



NF_EXT = 3 

LAMBDA_5 = 0.226232 GEV NLO-ALPHA_S (M_Z) = 0.118000 

REN-SCALE = 500.000 GEV 

T-MASS = 176.000 GEV B-MASS = 5.00000 GEV 

MSSM (2-LOOP): TG(BETA) = 1.50 



Z-MASS = 91.187 GEV W-MASS = 80.330 GEV 

NO MIXING SUSY-SCALE = 1000.0 GEV 

HIGGS = H 



M_H 




500.000 


GEV 






M_A 




490.788 


GEV 






CH.B 




1 . 47436 


G~H. 


.T 


-0.691622 


GAM.LO 




9.15160 


MEV 






E.VIRT 




15.2434 








E_GGG 




7.68668 


+- 


0, 


,812334E-04 


E_GQQ 




-0.972425 


+- 


0, 


.307239E-04 


E_TOT 




20.0128 


+- 


0, 


. 122860E-03 


GAM.NLO 




13.9932 


+- 


0, 


. 297226E-04 MEV 
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